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ABSTRACT 

We measure the relative evolution of the number of bright and faint (as faint as 0.05 
L*) red galaxies in a sample of 28 clusters, of which 16 are at 0.50 ^ z ^ 1.27, 
all observed through a pair of filters bracketing the 4000 A break rest-frame. The 
abundance of red galaxies, relative to bright ones, is constant over all the studied 
redshift range, < z < 1.3, and rules out a differential evolution between bright 
and faint red galaxies as large as claimed in some past works. Faint red galaxies are 
largely assembled and in place at z = 1.3 and their abundance does not depend on 
cluster mass, parametrized by velocity dispersion or X-ray luminosity. Our analysis, 
with respect to previous one, samples a wider redshift range, minimizes systematics 
and put a more attention to statistical issues, keeping at the same time a large number 
of clusters. 

Key words: Galaxies: evolution — galaxies: clusters: general — Galaxies: luminosity 
function, mass function — Galaxies: formation 



1 INTRODUCTION 

The evolution of faint red galaxies in clusters is a highly 
debated topic for two reasons: different observers have 
claimed controversial results, and clusters of galaxies are of- 
ten claimed to be interesting laboratories where studying 
the effect of the environment. Red galaxies, in particular, 
have different assembly histories in halos of different masses, 
yet observationally the detection of a environmental depen- 
dence of their properties escapes a detection. For example, 
differences between cluster and field fundamental planes are 
small, if any (Pahre et al. 1998), so small that the Coma clus- 
ter fundamental plane (Jorgensen et al. 1996) is routinely 
used as zero-redshift reference for studying the evolution of 
field galaxies, and so small that previously claimed differ- 
ences are probably due to having overlooked the difficulty 
of the statistical analysis (van Dokkum & van der Marel 
2007). Similarly, the colour of the red sequence seems not 
to depend on clustercentric distance (Pimbblet et al. 2002; 
Andreon 2003) or galaxy number density (Hogg et al. 2004, 
Cool et al. 2006). 

The red colour, by which red galaxies are defined and 
selected, induces a selection effect: at every redshift only 
galaxies whose stellar populations are red (i.e. old, mod- 
ulo dust, of no interest here) enter the sample. It is not a 
surprise, then, to find old-selected populations to be old. A 
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different question is whether galaxies that have an old stellar 
population were fully assembled at early or late times. An- 
swering this question requires a measurement of the abun- 
dance of red galaxies as a function of look-back time. For 
clusters, there is a further complication: clusters have differ- 
ent richnesses, jeopardizing any look-back time trend if the 
richness dependence is not factored out. It is easy, further- 
more, to qualitatively claim that the red sequence is built 
later (i.e. a lower redshift) in poor environments than it is 
in dense environments, but this might just be do to signal 
to noise issues, because in poorer environments the red pop- 
ulation is a minority one, and its contrast with respect to 
other populations (e.g. background) noisier. A sound sta- 
tistical assessment of the abundance of faint red galaxies is 
therefore compelling. 

Usually, the richness dependence of the abundance of 
faint red galaxies is removed by normalizing it to the num- 
ber of bright red galaxies, i.e. by computing the faint-to- 
luminous ratio, or any related quantity, like the faint-end 
slope a of the luminosity function. The analysis of the faint- 
to- luminous ratio, performed by Stott et al. (2007), or its 
reciprocal, the luminous-to-faint ratio by De Lucia et al. 
(2007), both suggest an evolution of the relative abundance 
of faint red galaxies, in the sense that at high redshift there 
is a deficit of faint red galaxies per unit bright galaxy. On 
the other end, Andreon (2006a) suggests no deficit of red 
galaxies, using a very small cluster sample, and Andreon et 
al. (2006) discard a considerable build up of the red sequence 
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Table 1. The ACS z 0.5 cluster and control field samples 



Name 


z 




Filters 








blue 


red 


T V WJ 
J-j V llA. VV 


1 27 


3 


F775W 


F850LP 


T .irnv T-T. 

i-J V li-A. Ill 


1 26 


3 


F775W 


F850LP 


RDCS J1252-2927 


1.23 


4 


F775W 


F850LP 


RDCS J0910+5422 


1.11 


1 


F775W 


F850LP 






]^ 


F606W 


F814W 




0.90 


]^ 


F606W 


F814W 


IWGA J1226. 9+3332 


0.89 


6 


F606W 


F814W 


MACS J0744.8+3927 


0.70 


1 


F555W 


F814W 


MACS J2129. 4-0741 


0.59 


1 


F555W 


F814W 


MACS J0717.5+3745 


0.55 


1 


F555W 


F814W 


MACS J1423. 8+2404 


0.54 


1 


F555W 


F814W 


MACS J1149. 5+2223 


0.54 


1 


F555W 


F814W 


MACS J0911. 2+1746 


0.50 


1 


F555W 


F814W 


MACS J2214.9-1359 


0.50 


1 


F555W 


F814W 


MACS J0257. 1-2325 


0.50 


1 


F555W 


F814W 


CT344 




1 


F606W 


F814W 


B0455 




1 


F555W 


F814W 


GOODS+PAN 




- 30 


F775W 


F850LP 


1 number of ACS field of 


view per 


filter. 







All clusters have coordinates and redshift listed in NED, except 
for MACS clusters, hsted in Ebehng et al. (2007). MACS clusters 
have been also studied by Stott ct al. (2007). 



on the basis of fossil evidence. Evidences presented in ear- 
lier works have been discussed in the mentioned papers and 
references therein. 

In this paper, we aim to understand if the colour- 
magnitude relation has been build up at early or late times, 
by studying many galaxy clusters at several look-back times. 

Throughout this paper we assume Q,m = 0.3, = 0.7 
and Ho — 70 km s~^ Mpc~^. All results of our stochastical 
computations are quoted in the form x ± a where x is the 
posterior mean and a is the posterior standard deviation. 



2 DATA & DATA REDUCTION 

This work makes use of archive data. The selection crite- 
ria used for the inclusion in our sample are the following: 
a) all observations must include a pair of filters bracketing 
the 4000 A break in the cluster rest-frame; b) control field 
observations with identical conditions (same telescope, in- 
strument, filters, dept and seeing) as cluster observations 
must be available; c) observations had to be deep enough 
to measure the faint end slope of the luminosity function; 
d) clusters had to be spectroscopically confirmed; e) data 
should be public available at the start of the work. 

Our sample is formed by three sets: a) 15 high red- 
shift clusters observed with the Wide Field Camera of the 
Advanced Camera for Surveys (hereafter ACS, Ford et al. 
1998, 2002) of Huhhle Space Telescope (HST, hereafter); b) 
two low redshift clusters observed by the Sloan Digital Sky 
Survey (SDSS); and c) two luminosity functions from litera- 
ture, one for a z ~ 0.25 cluster sample and one for one more 
high redshift cluster observed by HST but with the Wide 
Field Planetary Camera 2. 

Table 1 lists the ACS sample, formed by 15 clusters 
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Figure 1. Rest-frame A sampling of the adopted filters for the 
clusters studied in this work. The shaded (yellow) band marks 
the 4000 A break. All clusters have been observed in a pair of 
filters bracketing the 4000 A break. The sample is formed by 28 
clusters, some of which have very similar redshifts and do not 
show up individually in the Figure. 
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Figure 2. Minor role of non-Poisson variance. The left panel 
shows galaxy counts of red galaxies in two widely different sky 
direction, whereas the right panel shows their colour distribu- 
tion. Differences between sky directions arc within Poisson fluc- 
tuations. 



at 0.5 ^ 2: ^ 1.27. More than 150 HST orbits, devoted to 
clusters, have been reduced and analized for this paper. As 
we need to statistically discriminate against fore- and back- 
ground interlopers, Table 1 also lists the adopted control 
fields. A control field matching the filter pair used for clus- 
ters is available for all targets. 

In order to provide a local [z ~ 0) reference, we use 
SDSS u, g data of two nearby clusters: Abell 1656 (A1656 
hereafter, i.e. Coma, z = 0.023) and AbeU 2199 (A2199 here- 
after, 2 = 0.030). Given the large SDSS sky coverage, the 
control field for our nearby clusters is taken all around them. 

Finally, the luminosity function of ten z 0.25 clusters, 
observed in B and I by Small et al. (1998), and of MS 1054 at 
z ~ 0.8, observed with HST Wide Field Planetary Camera 
2 in F606W and F814W and presented in Andreon (2006) 
have been taken from the literature. These LFs are fully 
homogeneous to those computed in this work. 

Figure 1 shows that all clusters have a pair of filters 
bracketing the 4000 A break. 

The raw ACS data listed in Table 1 were processed 
through the standard CALACS pipeline (Hack 1999) at 
STScI. This includes overscan, bias, and dark subtraction, as 
well as flat-fielding. Image combination has been done with 
the multidrizzle software (Koekemoer et al. 2002). The data 
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Figure 3. Luminosity function of red galaxies in some (to save space) of our clusters. The solid curve and the shaded region mark the 
LF and its (highest posterior) 68 % error region as computed with Bayesian methods. Points and error bars mark simply-derived LFs, 
computed as difference of counts in the cluster and reference lines of sight, using "two u from the colour-magnitude relation" as definition 
of red. The top-left panel refers to the two Lynx clusters stacked together. 



quality arrays enable masking of known hot pixels and bad 
columns, while cosmic rays and other anomalies are rejected 
through the iterative drizzle/blot technique. Sources are de- 
tected using SExtractor (Bertin & Arnout 1996), making use 
of weight maps produced by Multidrizzle. Star/galaxy sep- 
aration is performed by using the stellarity index given by 
SExtractor. HST images are calibrated in the Vega system, 
using the zero points provided in the HST data handbook. 
Completeness is computed as in Garilli, Maccagni & An- 
dreon (1999), from the brightest luminosity of the detected 
objects of faintest surface brightness. Only data brighter 
than the completeness magnitude are kept. 

All science (i.e. cluster) and two of the control fields, 
CT344 and BO0455, have been combined (and catalog built) 
by ourself, while the remaining control field, GOODS+PAN, 
has been generously given to us by D. Macchetto. These im- 
ages come from the same telescope, instrument and filters 
and have been processed with the same software as science 
data (i.e. CALACS, Multidrizzle and SExtractor), but have 
been combined by someone else (than the author) . By reduc- 
ing by ourself part ot the GOODS-I-PAN data, we checked 
that their and our reductions are indistinguishable. 

The left panel of Figure 2 shows galaxy counts of red 
galaxies, where 'red' is taken to mimic our later selection, 
for two widely different sky directions: CT344 and a field 
in Benitez et al. (2003). The right panel shows the colour 
distribution in the two directions. Differences between the 
two sky directions are comparable to Poisson errors on the 
average value. Therefore, for areas, magnitudes and colours 
of interest in this paper, non-Poisson fluctuations of galaxy 
counts can be neglected. 

For the nearby cluster sample, catalogs have been ex- 
tracted from the SDSS 5th data release (Adelman-McCarthy 
et al., 2007), which have been produced by the SDSS pipeline 



and are not calibrated in the Vega system. We checked that 
synthetic f/, V and U — V computed from u, g SDSS photom- 
etry is indistinguishable from observed U, V and U — V pho- 
tometry for red galaxies in A1656 cluster direction, taken 
from Terlevich et al. (2001), and derived with traditional 
techniques (stare exposures, calibration in the Vega system, 
and catalogs built with SExtractor). 



3 FAINT-TO-LUMINOUS RATIO OF RED 
GALAXIES 

We modelled the distribution of galaxies in the red se- 
quence as Gauss-distributed in colour at every magnitude 
and Schechter (1976) distributed in magnitude. The mean 
colour of the Gauss varies linearly with magnitude, because 
the color-magnitude relation is linear. Furthermore, we allow 
a broadening of the colour-magnitude relation due to both 
photometric errors and an intrinsic scatter. As explained in 
appendix A in Andreon (2006a), with Bayesian methods we 
solved at once for all parameters (colour-magnitude slope, 
intercept and intrinsic scatter, characteristic magnitude M*, 
faint-end slope a, and normalization of the Schechter, back- 
ground parameters), hence fully accounting for the back- 
ground (including uncertainty, variance and covariance with 
all parameters). 

Our definition of red is "galaxies under the Gauss cen- 
tered on the red sequence", similar to some SDSS works 
(e.g. Balogh et al. 2004; Ball et al. 2006). Previous studies 
(Andreon et al. 2006, De Lucia et al. 2007) have shown that 
the precise definition of 'red' has a negligible impact on the 
results. We have checked it for our own sample and our def- 
inition, by adopting a simpler definition of red (within 2a 
from the red sequence, plotted as dots in Fig 3). 

The luminous-to-faint ratio is computed as the ratio 
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of the number of galaxies on the red sequence in appropri- 
ate absolute magnitude ranges. The number of galaxies in 
a given range is, by definition, the integral of the luminos- 
ity function over the concerned range. The range definitions 
are taken from De Lucia et al. (2007): Mv < —20 mag 
and —20 < My < —18.2 mag. Magnitudes are passively 
evolving, modelled with a simple stellar population of solar 
metallicity, Salpeter IMF, from Bruzual & Chariot (2003), 
as in De Lucia et al. (2007). As a sanity check, the same 
model has been checked to reproduce the colour of the red 
sequence at My ~ —20.0 mag for all our clusters. 

From now on, the two Lynx clusters are stacked to- 
gether to improve the signal to noise. LFs are computed for 
galaxies within a cluster-centric radius listed in Table 2. The 
considered region has been chosen as a compromise between 
sampling a large portion of the cluster and not including 
a too large contribution from background galaxies. MACS 
clusters are larger that the instrument field of view, and 
therefore we choose the largest radius that fit in the fully 
exposed part of the image, consistently with the choices of 
Andreon (2006) and Small et al. (1998), whose LF are in- 
cluded in the present work, as mentioned. 

The potential dependency of the LF slope on the con- 
sidered cluster portion has a small impact on our study, 
because we explicitely allows the observed value of the slope 
to scatter around to its true value by more than its uncer- 
tainty. In fact, in the next section we we allow an intrinsic 
scatter in our model: see eq. 1. This argument is developed 
further in Sec 4.2. 



4 RESULTS 

Figure 3 shows the luminosity function of a sub-set (to save 
space) of studied clusters. 

Figure 4 shows the luminous-to-faint ratio, L/F, as a 
function of redshift for clusters with LF measured in this 
paper (solid points). Our data are in agreement with De 
Lucia et al. (2007) data (open points), but our wider redshift 
coverage suggests a shallower trend than the one hinted in 
De Lucia et al. (2007) from their data points. 

In this work we refrain to perform inferences using L/F 
or its reciprocal, F/L, for reasons detailed in sec 4.4, mainly 
of statistical nature. The use of the faint end slope, a, is a 
measure of the faint-to-luminous ratio, it is easier to deal 
with from a statistical point of view, and has the further 
advantage that it uses all the data available, including data 
fainter than —18.2 mag that would be otherwise wasted us- 
ing L/F. 

Figure 5 shows the slope, a, as a function of redshift 
for the whole cluster sample, i.e. for 28 clusters, of which 16 
at z ^ 0.5. Marginalization accounts for the known correla- 
tion between parameters (e.g. M* and a). For example, the 
large error of some data points is due to the fact that many 
(M*,a) pairs fit almost equally well the data and thus a large 
range of a values is acceptable, a errors also account for dif- 
ferences in the galaxy background counts in the cluster and 
control field lines of sight, because, as mentioned, we "solve" 
for all parameters at once (technically, we marginalize over 
other parameters). Table 2 lists the a and L/F values found. 

The data are in agreement with the lack of a deficit of 
faint red galaxies suggested by Andreon (2006) on the basis 
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Figure 4. Relative abundance of faint and bright red galaxies, 
as parametrized by the luminous-to-faint ratio, for clusters with 
LF measured in this paper (solid dots) and in De Lucia et al. 
(2007, open points). Although in agreement, our data indicate a 
shallower trend with redshift than De Lucia et al. (2007) data 
points. Two points at z = 0.55 fall one on the top of the other. 
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Figure 5. Relative abundance of faint and bright red galaxies, as 
parametrized by the faint slope a of the cluster luminosity func- 
tion for the whole sample of 28 clusters studied here. The red- 
shift dependence of the relative abundance of bright red galaxies 
is small, if any. The point at z = 0.25 (z = 1.27) is the average 
of 10 (2) clusters. The shaded (yellow) region shows the (highest 
posterior density) 68 % error region. The dashed lines delimit the 
itlo" intrinsic (i.e not accounted for measurement error) scatter. 
The trends proposed by Stott et al. (2007) and we fitted on De 
Lucia et al. (2007) data are also marked with solid lines (labelled 
by "S07" and "DL07", respectively) up the the largest studied 
redshift by them, and marked with dotted lines afterward. 



of a very small sample of clusters and reject some trends sug- 
gested in previous works. Let's consider: a) our maximum 
likelihood fit of the L/F data points in Fig. 9 of De Lucia 
et al. (2007) and b) the Stott et al. (2007) F/L vs z fit. 
The two fits have been transformed in q vs z trends using 
the L/F, F/L and a definitions. Figure 5 shows that at low 
and intermediate redshift the De Lucia et al. (2007) trend, 
marked with "DL07" , is compatible with our data. However, 
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a constant, i.e. a more economical model having one degree 
of freedom less, also well describes our data (and also theirs, 
see Fig. 5) over the common redshift range {z < 0.8) and, 
actually, also above. Furthermore, neither De Lucia et al. 
(2007) nor our data request a more complex model than a 
constant plus an intrinsic scatter. The computation of the 
Bayes factor shows that the De Lucia et al. (2007) trend 
is disfavored, with respect to 'no trend at all' by our data 
with odds 14:1, i.e. there is moderate evidence against an in- 
crease of the luminous-to-faint ratio as large as pointed out 
by De Lucia et al. (2007). We refrain, therefore, from fitting 
a more complex model, and we adopt a constant model. Fig- 
ure 5 also shows that the Stott et al. (2007) fit, marked with 
"S07" in the Figure, nicely reproduces the observed values 
in the reduced redshift range, 0.5 ^ 2; ^ 0.6, where we share 
clusters and HST data with them, but disagrees outside it, 
in particular at low redshifts. Furthermore, in the local uni- 
verse, the Stott et al. (2007) fit and data also disagrees with 
De Lucia et al. (2007) data and trend. Our data clearly dis- 
card the trend proposed by Stott et al. (2007). 

Using Bayesian methods (D'Agostini 2003, 2005) and 
uniform priors we 'fitted' the data point with a constant, 
accounting for errors and allowing an intrinsic (i.e. not ac- 
counted for errors) Gaussian scatter, A/'(0, (Jintr)- We found: 
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y{z) = -0.91 ± 0.06 + A/'(0, 0.13 ± 0.06) 



(1) 



displayed in Fig 5. Using our own data alone, i.e. ignor- 
ing the Small et al. (1998) z = 0.25 composite cluster, we 
found an indistinguishable result: 

a{z) = -0.89 ± 0.06 + A/'(0, 0.16 ± 0.06) (2) 



4.1 Richness dependency 

Koyama et al. (2007), on the basis of three clusters at 
z ^ 0.8, suggested that the relative abundance of faint red 
galaxies is dependent on cluster richness or mass (actually. 
X-ray luminosity in their work), in the sense that poorer 
systems show stronger deficits. Fig. 6 and 7 plot two deficit 
estimators, the slope a and the L/ F ratio, vs two mass es- 
timators. X-ray luminosity and velocity dispersion, for clus- 
ters at z ^ 0.5. X-ray luminosities are taken from Ettori 
et al. (2004), Lubin et al. (2004) and Ebeling et al. (2007). 
They comes from Chandra or XMM pointed observations 
and have, typically, errors of 10 % or less. MACS luminosi- 
ties, in the 0.1-2.4 keV band are converted in bolometric 
assuming a thermal bremsstrahlung spectrum with the mea- 
sured temperature. Velocity dispersions are taken from Ebel- 
ing et al. (2007), Stanford et al. (2001) and Maughan et al. 
(2004). Solid dots emphasize clusters in the reduced redshift 
range 0.50 ^ z 0.70, to limit the (negligible, see previ- 
ous section) efi'ect of evolution. There is no obvious trend 
between cluster mass and the relative abundance of faint 
red galaxies. (Green) crosses, we connected by a solid line 
in the top panel of Fig 6, show the three clusters studied 
by Koyama et al. (2007), one of which is MS1054.4-0321, 
the latter taken from Andreon (2006). The slope a is de- 
rived by us from their luminous over faint ratio assuming 
a Schechter function. The figure shows that the observed 
slopes a are plausible, since two other clusters in our sam- 
ple show similar values of the relative abundance of faint red 
galaxies. However, the trend suggested by these three points 
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Figure 6. Cluster mass dependency, as parametrized by X-ray 
bolometric luminosity, of the relative abundance of faint red 
galaxies, as parametrized by the a (top panel) or the luminous- 
to-faint ratio (bottom panel), for clusters with z ^ 0.5. Clusters 
in the 0.5 ^ ^ ^ 0.70 range are marked with a solid (red) dot. The 
three crosses connected by a solid line in top panel mark clusters 
used by Koyama et al. (2007) to suggest a mass-dependent trend 
in the relative abundance of faint red galaxies. 



(slanted line) is clearly too steep, and obviously ruled out 
by our data. Finally, we note that RXJ0152. 7-1357 point, 
i.e. the middle point of the three plotted, has been put in 
Koyama et al. (2007) and in our plot at the sum luminos- 
ity of the two sub-clumps that form the cluster, not at the 
mean L/F value. Would the RXJ0152. 7-1357 point be put 
at the average x-ray luminosity of the two clumps, the one 
typically experienced by galaxies in this clusters and con- 
sistently with the choice of quoting a mean L/F, the three 
L/F points would no longer show any monotonic trend with 
X-ray luminosity. 

In conclusion, the abundance of faint red galaxies does 
not considerably depend on cluster mass (in the range sam- 
pled by data, of course) with the Koyama et al. (2007) trend 
largely based on a sample of inadequate size, given the large 
intrinsic scatter, a possibility also mentioned by these au- 
thors. 
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Figure 7. As Figure 6, but using velocity dispersion as cluster 
mass estimator. 

4.2 Radius- and scatter- related effects and the 
advantages of allowing an intrinsic scatter 

Clusters have no sharp boundaries. In order to understand 
the potential effect of the choice of the studied cluster por- 
tion, we compute ct and LjF oi A1656 cluster within two 
clustercentric radii: 0.2 and 0.8 deg. A1656 cluster has been 
chosen because it has the best determined values of a and 
LjF among all our clusters. For A1656, the relative abun- 
dance of faint red galaxies, as parametrized by these quan- 
tities, is the same within the two radii. Although this test is 
reassuring, we cannot generalize from a single example. 

Our model, eq. 1, explicitly allows an intrinsic variance 
in the relative abundance of faint red galaxies, due for ex- 
ample to the mentioned clustercentric dependence, cluster- 
to-cluster or other possible (unknown, for the time being) 
systematics. These terms are thus a source of scatter, not 
different from a random number added to each measure- 
ment. We model such random process, whatever its phys- 
ical nature is related to clustercentric distance, to cluster- 
to-cluster variance or whatsoever unidentified reason, with 
a normal (Gaussian) of unknown variance (eq. 1), for lack 
of evidence toward any more complex model. In passing, 
in absence of more information, the Gaussian is the maxi- 
mum entropy choice among all real- valued distributions with 
specified mean and standard deviation (e.g. Sivia 2006). The 
sum rule of probability states that in order to proceed with 
the inference, we need to marginalize over this (nuisance) 
parameter. The intrinsic scatter parameter, and its conse- 
quent marginalization, offers protection against claiming a 



trend when just too few data are available: marginalization 
spreads the probability of a (redshift or mass) trend over 
a large range of slope values, i.e. quantifies the researcher 
good sense that when just a few data points are available 
and an intrinsic scatter is there, one should be prudent in 
claiming the existence of a trend. Assuming a single value 
for the intrinsic scatter, as other authors sometime implic- 
itly take when looking for a trend, artificially collapses the 
error ellipse along this axis and leads to determinations with 
overly-optimistic confidence. 

Therefore, since the model allows an intrinsic scatter, 
the analysis of the redshift depencence of the relative aboun- 
dance of faint red galaxies is correct even if: a) a cluster- 
centric dependence of the the relative abundance of faint 
red galaxies exists, provided that we are not sampling in- 
creasing smaller cluster portions as the redshift increases or 
decreases, or b) the relative abundance of faint red galaxies 
differs from cluster to cluster. Both cases are just source of 
scatter and our model account for them. 

As mentioned, we found a non-zero intrinsic (i.e. not 
accounted for measurement errors) scatter [uintr = 0.13 ± 

0. 06), quantifying past claims of a heterogeneity in the rela- 
tive abundance of faint red galaxies. Inspection of the colour- 
magnitude relations of outliers in Fig 5, e.g. MACS J0257.1- 
2325, confirms that these clusters have an underpopulated 
red sequence at faint magnitude or an overpopulated one at 
bright ones. 

Therefore, in the quest of a build-up of the red sequence, 
an intrinsic scatter must be allowed, in order not to over- 
weigh 'outlier' clusters, and not to overstate the precision 
and the statistical significance of the found (redshift, mass 
or whatever) trend. We note that the existence of an intrin- 
sic scatter has been claimed in previous works (De Lucia et 
al. 2007, Stott et al. 2007) but ignored (Stott et al. 2007) 
or not rigorously accounted for (De Lucia et al. 2007) when 
establishing the veracity of the claimed redshift trend. 

The existence of an intrinsic scatter testifies that: a) 
there is a yet to be identified physical mechanism that affects 
the relative abundance of faint red galaxies, and b) current 
data are of adequate quality to perform such measurement, 

1. e. that the topic deserves further investigation. 

4.3 Joining high- and low- redshift information 

We emphasize that our (past) knowledge about individual 
nearby clusters tell us that at least some galaxies on the 
red sequence have a spiral morphology (Butcher & Oemler 
1984; Oemler 1992; see Fig. 3a in Andreon et al. 1996 or Fig. 
4 in Terlevich et al. 2001 for A1656 galaxies). Their spiral 
arms testify that, in the past, these galaxies were forming 
stars, i.e. were blue, and therefore were not on the red se- 
quence. Furthermore, at least for A1656 (Coma), red spirals 
have lower surface brightness than blue spirals (Andreon 
1996), as expected if the former are the descendent of the 
latter. Since, on average, spirals are fainter than early-type 
galaxies (e.g. Bingelli, Sandage, Tammann 1988 for Virgo, 
Andreon 1996 for A1656), we expect that the abundance 
of faint red galaxies grows somewhat with time, just be- 
cause of the evolving colour (toward the red) of some spi- 
rals. However, it cannot grow too much, otherwise it would 
bend the colour-magnitude and inflate its scatter. The ar- 
gument is the usual one (e.g. Bower, Lucey & Ellis 1992): 
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Table 2. Extraction radius r in arcmin, faint end slopes a and 
luminous-to- faint ratios L/F 



Cluster name 


r 


a 


L/F 


Lynx E+W 


1.0 


-1.12 + 0.22 


0.57 + 0.15 


RDCS J1252-2927 


1.0 


-0.89 + 0.37 


0.57 + 0.16 


RDCS J0910+5422 


1.0 


-0.77 + 0.36 


0.60 + 0.14 


GHO 1602+4329 


1.0 


-0.54 + 0.38 


1.03 + 0.42 


GHO 1602+4312 


1.0 


-0.67 + 0.33 


0.84 + 0.23 


IWGA J1226. 9+3332 (CL1226) 


1.5 


-0.97 + 0.33 


0.29 + 0.07 


MACS J0744.8+3927 


1.4 


-0.34 + 0.31 


0.67 + 0.12 


MACS J2129.4^0741 


1.4 


-0.78 + 0.18 


0.52 + 0.08 


MACS J0717.5+3745 


1.4 


-1.04 + 0.07 


0.72 + 0.06 


MACS J1423.8+2404 


1.3 


-0.89 + 0.16 


0.74 + 0.14 


MACS J1149. 5+2223 


1.2 


-0.92 + 0.09 


0.65 + 0.07 


MACS J0911. 2+1746 


1.2 


-0.78 + 0.16 


0.53 + 0.08 


MACS J2214.9-1359 


1.3 


-0.76 + 0.12 


0.71 + 0.09 


MACS J0257.1-2325 


1.3 


+0.02 + 0.31 


0.35 + 0.07 


A2199 


48 


-1.12 + 0.06 


0.48 + 0.04 


A1656 


48 


-1.02 + 0.05 


0.44 + 0.04 



a) a heterogeneity in the star formation history leads to a 
heterogenous population in colour (unless something else co- 
ordinately conspires to keep the colour scatter small) ; and b) 
a delayed stop of the last star formation episode delays the 
arrival of a galaxy on the colour-magnitude relation, bending 
it (or increasing the colour scatter if there is an un-delayed 
population). In Abell 1185 the colour-magnitude relation is 
linear and the scatter in colour is small (0.04 mag) down to 
M* + 8 (Andreon et al. 2006). In A1656, the scatter stays 
con.stant to low levels (0.0-5 mag) down to M* + 4 (Eiscn- 
stein et al. 2007). Therefore, fossil evidence points towaxd 
a small, but not null, differential build up of red sequence 
galaxies. 

Theory (De Lucia et al. 2007) shows that a model in 
which star formation histories of blue galaxies are truncated 
produces an important change in the luminous-to-faint ratio, 
larger than allowed by De Lucia et al. (2007) data, producing 
too many faint galaxies by a factor two (but note that these 
authors consider it as 'approximatively consistent' with their 
data). Direct measurement of the abundance of red faint 
galaxies over redshift (this paper) indicates a shallower trend 
with redshift than suggested by theory at a point that data 
are consistent with no trend at all. 

Therefore, data at cosmological redshifts and the tight- 
ness of the colour-magnitude relation at low redshift strongly 
argue against a scenario where many blue galaxies trans- 
form themselves in faint red galaxies, whereas the presence 
of some spiral galaxies on the red sequence in nearby clusters 
suggests a redshift trend in the relative abundance of faint 
red galaxies should be observed in sufficiently large cluster 
samples, although none is clearly revealed in the present one. 



4.4 Some advantages and shortcomings of current 
studies 

Other authors argue that L/F or its reciprocal, F/L, are 
preferable to a in the study of the relative abundance of faint 
red galaxies, usually with the rationale that the Schechter 
function might not describe the luminosity function in the 
studied mag range. Beside the fact that the very same au- 



thors find acceptable fit on their data (typically, x2 ~ 1 
values), and thus they argue something not supported by 
their own data, we note that L/F, and its reciprocal, F/L, 
both are quantities difficult to manage from a statistical 

■^oint of view. For example an average value, computed by 
a weighted sum, or a fit performed minimizing the x^i has 
a special meaning, because the result depends on whether 
L/F or F/L is averaged (fitted). For example, let's consider 
two, for sake of clarity, data points, {f /l)i = 3 ± 0.9 and 
(//O2 ~ 0.3.33.3 ± 0.1 and two possible averages. The error 
weighted average {f /I) is 0.37. The error weighted average 
(l/f) of the reciprocal values {{l/f)i = !/(///),; 0.3333±0.1 
and 3 ±0.9) is again 0.37, fairly different from the reciprocal 
of (f/l), 1/0.37 = 2.7. Therefore (///) / 1/(1/ f)- At first 
sight, by choosing the parametrization of the aimed quantity, 
the astronomer may chose the result he want. Furthermore, 
(f/l) has a value near to the point with index 2, (//O2, 
whereas {l/f} has a value near to the other data point, with 
index 1, {f/l)i, a strange situation, indeed. Similar prob- 

Jems are present with two data points diflFering by just la, 
or with small samples. Therefore, astronomers who want to 
use L/F or F/L are invited first to understand what is go- 
ing on in the simple case of just two measurements and a 
fit with a constant (the mentioned weighted average), and 
then proceed to the case they are really interested in: a few 
data points and a fit with one more degree of freedom (the 
redshift or mass dependence). 

As mentioned, there are currently a few determinations 
of the evolution of faint red galaxies. Some of them make 
different claims concerning the deficit of faint red galaxies at 
high redshift, yet we have verified that often, but not always, 
data agree with each other, as shown in Fig. 4 and Fig. 5, 
in the (usually small) range where they are well determined. 
The present work offers some advantages with respect to 
previous ones. 

- First, our determination is more sensitive to evolu- 
tion, because our cluster sample displays the widest redshift 
coverage while keeping a large sample of clusters. 

- Second, the present work minimizes systematics, for 
example using a colour index bracketing the 4000 A break at 
every redshift. Fig 1 shows that SDSS u and g filters and, at 
a lower extent, B and I filters aX z = 0.25, sample the 4000 
A break in similar way as HST filters do at higher redshift. 
This is not the rule: for example, De Lucia et al. (2007) use 
rest-frame U — V at high redshift but r' — i' at z = for 
C4 clusters (r' and i' have effective A = 6165, 7481 A, re- 
spectively). The central wavelength of the bluest filter (R) 
used in the recent 0.35 < 2 < 0.95 study by Gilbank et al. 
(2007) lies longward of the Balmer break for z < 0.75 (Fig 
8). By using the same filter pair {R — z') at all redshifts, 
rather than a color selection mimicking U — V at sl\ red- 
shifts, their technique introduces a potential bias, as several 
spiral types move from red to blue between the low and the 
high redshift samples because of their observational strat- 
egy. Furthermore, at 2; > 0.55 the data used by Gilbank 
et al. (2007) are not deep enough to sample faint galaxies 
{My = —18.2 mag, such as those considered here, in Stott 
et al. 2007, de Lucia et al. 2007, in Gilbank & Balogh 2008, 
etc.), see top panel of Fig 8. In their later work, Gilbank 
& Balogh (2008) use data lacking appropriate Balmer break 
coverage (the z = 0.25 point from Hansen et al. 2008 and 
the 2 = point from De Lucia et al. 2007) and omit data 
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Figure 8. Bottom panel: Rest-frame A sampling of the filters used 
in Gilbank et al. (2007). The shaded (yellow) band marks the 
4000 A break. The central wavelength of the bluest filter (R) lies 
longward of the Balmer break for z < 0.75. Compare with Fig 1. 
Upper panel: R and z' mag of a passive evolving My = —18.2 mag 
galaxy (red dashed and blue solid curves, respectively) and the 
limiting magnitude of the observational material used by Gilbank 
et al. (2007). Gilbank et al. (2007) data arc not deep enough to 
detect these faint galaxies at 2 > 0.55. 



with better 4000 A sampling (Fig 1) and appropriate dept 
(those in the present work). In passing, we also disagree with 
their statement about the number of clusters at low redshift 
in our sample: their claim that it "only contains two z < 0.5 
clusters", while instead our sample includes 12 clusters, 10 
from Small et al. (1998) and two from our own analysis of 
SDSS data. 

- Third, interlopers are removed using observations 
taken in the very same bands as cluster observations (see 
Table 1), to avoid systematics (see Small et al. 1998 and 
Andreon et al. 2006). This is often not the case: for clus- 
ter and control field Stott et al. (2007) use different filters, 
whereas de Lucia et al. (2007) use different telescopes and 
filters. The impact of these systematics is not quantified in 
the mentioned works. 

-Fourth, we feel our statistical analysis to be preferable: 
beside already mentioned statistical considerations, there 
are a number of debatable issues in other works, such as 
averages of incompatible measurements, Poisson errors for 
binomial distributed quantities, and unphysical results such 
as negative number of galaxies. 

- Finally, our clusters are spectroscopically confirmed 
and have an X-ray emission that confirms the existence of 
deep potential wells. Instead, we ignore whether candidate, 
or putative, clusters without a spectroscopic confirmation or 
an X-ray detection, studied in some other papers (e.g. Ko- 



dama et al. 2004) are clusters or line of sight superpositions. 
Sometimes, foUowup spectroscopic observations shows that 
a considerable fraction of them are line of sight superposi- 
tions (e.g. Yamada et al. 2005). 



5 SUMMARY AND CONCLUSIONS 

The history of mass assembly of bright (massive) red galax- 
ies in clusters is pretty well known: they were assembled 
at early times, as testified by the passive evolution of their 
characteristic magnitude (e.g. de Propris et al. 1998, 2007; 
Andreon et al. 2007), the constancy of their stellar mass 
function (e.g. Andreon 2006b) and of the halo occupation 
number, i.e. the number of galaxies per unit cluster mass 
(Lin et al. 2006, Andreon et al. 2008). We stress that aU 
mentioned works favour the above scenario, but only one, 
(Andreon 2006b), excludes contender models, and we em- 
phasize that most mentioned works have samples that are 
dominated, but not exclusively composed, by red galaxies. 

The history of mass assembly of faint red galaxies is far 
less clear. The present paper shows that a non-evolution of 
the faint end slope a, or any related number such as the 
luminous-to-faint ratio, is fully compatible with the data. 
This implies that the history of mass assembly of faint 
red galaxies is strictly parallel to the one of their massive 
cousins, in order to keep the relative abundance constant. 
Therefore, the build up of the red sequence is largely com- 
plete by z = 1.3 down to 0.05L* , and, if a differential filling 
is envisaged, it should occur mostly at much larger redshift. 
Similarly, cluster mass, as parametrized by X-ray luminosity 
or velocity dispersion, seems not to play any role in shaping 
the relative abundance of faint galaxies, contrary to some 
previous claims. Our claims are based on one of the largest 
samples, spread over the wider redshift range studied thus 
far with a large cluster sample, with great attention to sys- 
tematics. A recent {z < 1.3) transformation of many blue 
galaxies in faint red galaxies would modify the faint-end 
slope of the LF, change the F/L ratio and infiate the color 
scatter of the colour-magnitude relation, none of which have 
been observed. Yet, a redshift trend is expected because of 
the spiral morphology of some faint red galaxies in nearby 
clusters, but a larger sample of clusters (at z ^ 0.5) is 
needed to measure its small amplitude. The present sam- 
ple is however large enough to discard the claimed steep 
trends previously suggested in literature. 
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